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Abstract 

We report the magnetic dipole and electric quadrupole hyperfine coupling constants of 6P3/2 state 
of ^^^Cs(I=|) obtained from the relativistic coupled cluster (RCC) method. To our knowledge, no 
prior electric quadrupole hyperfine CC calculation is available for ^ Cs(6P3/2) state. Our computed 
r^ ■ magnetic dipole (A) and electric quadrupole {B) coupling constants are in excellent agreement with 

^ 1' the experiment. 

PACS number(s). : 31.15.Ar, 31.15.Dv, 31.25.Jf, 32.10.Fn 

Theoretical determination of hyperfine coupling constant is one of most stringent tests of atomic wave 
functions near the nucleus. Also accurate prediction of hyperfine coupling constants and other related 
^ I properties (viz. transition energies) requires precise incorporation of relativistic and higher correlation 

c/3 ■ and relaxation effects as these effects are strongly entangled. In particular, the hyperfine coupling 

constant calculations are relevant to the parity non-conservation study in atoms because the electro- 
£^ ■ weak interaction is also a short range force. 

^H ^l The relativistic and dynamical electron correlation effects can be incorporated in many-electron sys- 

O ■ tems through a variety of many-body methods. The coupled cluster method (CCM) has emerged as 

one of the most powerful and effective tool for a high precision description of electron correlations in 

many-electron systems. The CCM is an all-order non-perturbative scheme, and therefore, the higher 

' ^ I order electron correlation effects can be incorporated more efficiently than using the order-by-order dia- 

(yr^ ■ grammatic many-body perturbation theory (MBPT). The CC method is size-extensive, a property which 

00 I has been found to be crucial for an accurate determination of state energies, bond cleavage energies and 

^^ ■ related spectroscopic constants. Since the order-by-order MBPT expansion terms are directly related 

to the CC equations (as the latter is an all-order version of the former scheme), the CC results can be 

^p ■ improved by adding the important omitted diagrams with the aid of low order MBPT. 

f^ I Precise measurements of hyperfine splittings in 6P3/2 state of ^^"^Cs have been carried out by Tanner 

"^ ■ and Weiman p]. The magnetic dipole and electric quadrupole hyperfine coupHng constant of 6P3/2 

O I state have also been determined from this high precision experiment. In this article, we report the 

' tj^ • magnetic dipole (A) and electric quadrupole (B) hyperfine coupling constant of 6P3/2 state of Cesium 

^~>| atom obtained from the RCC method to access the performance and accuracy of the RCC scheme. We 

i^ ■ reiterate that accurate computation of hyperfine coupling constants requires precise incorporation of 

electron correlations and accurate description of the wave functions near the nucleus. 

The relativistic hyperfine energy W for states J — L ± ^ is given by 

><■ W-A(I J) I p 6(/-J)^ + 3(/.J)-2((/./)(J.J) 
H' W-A{I.J)+B 2/(/-l)2J(J-l) ' ^^> 



where / is the nuclear spin. A and B are, respectively, the coefficients of magnetic dipole and electric 
quadrupole contribution to the hyperfine structure. The radial matrix elements related to A and B are 
available elsewhere j2]. [Note that B=0 for J = ^.] 

We employ the straight forward extension of non-relativistic open-shell coupled cluster theory to the 
relativistic regime by adopting the no- virtual-pair approximation (NVPA) along with appropriate mod- 
ification of orbital form and potential terms .SJ. We begin with Dirac-Coulomb Hamiltonian (Ho) 
expressed in normal order 

Hd = Hn- {0\H\Q) = Y.{i\f\j) {ala,} + J ^ {zj\\kl) {ala]aiak} , (2) 

ij i,j,k,l 



where 

{ij\\kl) = {tj\ — \kl}-{tj\ — \lk}. (3) 

ri2 ri2 

The valence universal Fock space open-shell coupled cluster method is employed which begins with the 
decomposition of the full many-electron Hilbert space of dimension N into into a reference space Mo of 
dimension M ^ N, defined by the projector P, and its orthogonal complement Mq associated with the 
projector Q = 1 — P. A valence universal wave operator SI is then introduced which satisfies 

I*,) = r!|vl/f )), z = l,...,M, (4) 

where \"^\ ) and \'^i) are the unperturbed and the exact wave functions of the ith eigenstate of the 
Hamiltonian, respectively. The wave operator fl, which formally represents the mapping of the reference 
space A^o onto the target space Ai spanned by the M eigenstates |^i), has the properties 

np = n, pn = p, n^ = n. (5) 

With the aid of the wave operator fl, the Schrodinger equation for the M eigenstates of the Hamiltonian 
correlating with the Af-dimensional reference space, i.e. 

i/|*,)=-B,|^,), z = l,...,M, (6) 

is transformed into a generalized Bloch equation, 

Hnp = finrip = fiPH^^p, (7) 

where HqQ = PHflP is the effective Hamiltonian. Once Eq. IJSJ is solved for the wave operator fl, 
the energies Ei, i — 1, . . . , Af , are computed by diagonalizing the effective Hamiltonian -ffgff in the Af- 
dimensional reference space Mo- Following Lindgren's formulation of open-shell CC [^, we express the 
valence universal wave operator SI as 

n = {exp(5)}, (8) 

where S is the excitation operator and curly brackets denote the normal ordering. 

Once ^i is known, the one-electron properties (O) are evaluated from the following expression 

iO)f. - , ^^^'^'^-^ . (9) 

In the actual calculations, the ground and excited state properties of Cs are computed using the finite 
basis set expansion method (FBSE) [51 with a large basis set of Gaussian functions of the form 

f.,fc(r)=r'=-e-"-'-', (10) 

with fc = 0, 1, . . . for s, p, ... type functions, respectively. For the exponents, the even tempering condition 

a, = ao/3*-' (11) 

is applied. Here, TV is the number of basis functions for a specific symmetry. The starting point of this 
calculation is the generation of DF orbitals for Cs''" which defines the (0,0) valence sector of the Fock 
space. The ion is then correlated by CCSD and one-electron is then added, following the Fock-space 
scheme [S]: 

Cs+(0,0) + e — >Cs(0,l) 

Table 1 compares the magnetic dipole and electric quadrupole hyperfine coupling matrix element of 
^^'^Cs(I=7/2) obtained from the CC method with the experiment Q] and with linearized CC calculations 
[3 El- We also report the 6s -^ ^Vzj2 transition dipole matrix to demonstrate the consistency of this 
method. Our computed magnetic dipole and electric quadrupole hyperfine coupling matrix elements are 
quite accurate [off by 0.143 MHz for the magnetic dipole and 0.02 MHz for electric quadrupole case]. The 
CC transition dipole matrix element reported here is also in excellent agreement with the experiment 
[2]. Finally, we emphasize that to our knowledge no prior electric quadrupole hyperfine CC calculations 
are available for 6P3/2 state of Cs atom. 
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Table 1: Comparison of hyperfine {A and B) coupling (in MHz) and 6s^6p3/2 transition dipole matrix 
element of ^'^^Cs atom with the experiment and other correlated calculations. 



This work Others Experiment 

A{6P3/2) 50.419 49.7850 48.51|ai 50.275(3) T 

B(6P3/2) -0.559 -0.53(2) [ij 

65*1/2 -^ 6F3/2 6.372 6.3707 6.304|H| 6.36 [S| 
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